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1 Introduction 

Consider a proton at rest. Surrounding this proton there is a cloud of mesons 
{tt,p,u!,(J), f2, etc.), which is fairly diluted at a distance that is large compared 
to the proton radius. That is, every meson is well separated from the rest of the 
mesons. Furthermore, due to the low density of mesons, nuclear models of proton- 
meson interactions should work on this regime, making perturbation theory a valid 
approximation, because the expansions are not only on the powers of the proton- 
meson effective coupling, but the series is also suppressed by powers of the product 
between the meson mass and the physical meson-proton distance[|l|. 

Now let us imagine a high Q^ virtual photon in deep inelastic scattering, which 
due to the high Q^ has vanishing size. Sometimes this photon will collide with the 
proton core, which constitutes a typical deep inelastic scattering events, but in other 
instances it scatters off one of the isolated mesons in the disperse meson cloud. The 
photon then breaks down the meson, and the pieces of the broken meson fragment 
independently of what happens to the proton. In fact, if the meson is far away, and 
if for instance the meson is a neutral pion, the most probable outcome is that the 
proton core will not be affected by the interaction between the photon and the meson, 
and after the meson is broken, the proton core will maintain its identity as a proton. 
Of course, if on the other hand the meson is charged, or if the core suffers the effects 
from the hard interaction, the proton can get excited into final states such as neutron, 
a A(-,0,+,++), etc. 

The basic picture is then the following: after the interaction we have a broken 
meson and a baryon which basically does not move much. Now, let us boost this 
picture to the laboratory system of HERA. What we see is that the proton looses 
very little momentum and continues down the beampipe, and the meson fragments 
are observed in the central rapidity zone, with a rapidity gap between the meson 
fragments and the proton direction (forward direction of beampipe). 

This picture we just presented gives a natural simple explanation of the rapidity 



gap events observed at HERA0. It is the purpose of this paper to show that this is 
indeed the case, doing the corresponding calculation in detail. 

The peripheral meson model we are considering here is different from the Pomeron 
exchange model0], which is the one that is usually ascribed to these large rapidity 
gap events. In fact, these two models give quite different predictions. Specifically 
the meson model predicts the presence of several excited baryonic states, besides the 
proton, in the forward direction of HERA; while the Pomeron model predicts that 
only the proton will be present there. Therefore future measurements in these region 
should be able to clearly distinguish between the two models. 

2 The Model 

The model we are considering is based on the diagram shown in Fig. 1, where a 
rapidity gap appears between the baryonic beampipe system A^* and the hadrons 
comprising system X. 

The measured diffractive structure function is in general a function of three vari- 
ables, which can be taken as the photon momentum transfer Q^, the Bjorken variable 
Q'^/2p-q, and the variable j3 = Q'^/2q- {p—p')- Since the cross section is dominated by 
small t = {p — p'Y values (|t| ^ Q'^,Mx), both [3 and the variable x = q- {p — p')/p- q 
can be reconstructed from the measured quantities xbj = xf3 and Mx/Q"^ = {l—j3)/ (3. 

Therefore the whole process depicted in Fig. 1 implies that the diffractive struc- 
ture function can be expressed as: 

m 

Here fm/p{x, Q^) is the probability of finding the emitted meson carrying a fraction x 
of the proton momentum, and F^{j3, Q^) is the usual deep inelastic structure function 
of the corresponding meson. Then the photon scatters on a quark that carries a 
fraction j3 of the meson momentum. Although in principle the sum is over all mesons. 
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Figure 1: Diagram for the process e{k) +p{p) — > e(fc') +N*{p') +X, in which a meson 
produces a rapidity gap between A^* and X. 



we will see below that a limited number actually contributes to the process. Notice 
that we are in fact using an equivalent particle approximation, in which the process 
is dominated by the region in which the particle (in this case the meson) is close to 
its mass shell. 

Since we are not looking into spin effects, for simplicity we will assume that the 
partonic distributions inside a meson are not very different form those inside a pion^. 
Therefore our task is to work out expressions for the different distribution functions 
of mesons inside a proton. 

The equivalent meson approximation allows for a separation of the amplitude for 
the whole process into an amplitude for a transition to a baryon A^* and a near mass- 
shell meson, followed by the amplitude for the interaction of the meson with a particle 
a (in our case the off-shell photon). Using this it is easy to show that the external 



and internal cross sections are related through: 

dUap = fm/p{x, k^) dx daan, (2) 

where the meson distribution inside a proton is then: 

;„,^(.,,^.). /_^_Mf^,A. (3) 

1-x 

Here | M^p P is the squared unpolarized amplitude for meson- A^* emission (summed 
over final spins of the meson, and averaged over initial spins of the proton), t = —K^ 
is the momentum squared of the pion, and rrim and Mp are the meson and proton 
masses respectively. K^ is a momentum cut-off scale inherent to the formalism of 
equivalent particle approximation, and in our case it represents a parameter that 
fixes the separation of the pion from the color field of the proton. So we expect that 
its value be of the order of Aqcd or less, although in our calculation we will leave it 
as an adjustable parameter. Notice that, as expected, the integral is dominated by 
small t values. 

The amplitude Mmp contains a baryon-meson-baryon form factor, which can be 
taken as0]: 

where A is a parameter and mm is the meson mass. In the parameterization of Ref. 
[5], the exponent n is equal to 1, except for the pNA case, where n = 2.[^ 

2.1 Pion distribution inside the proton: 

We start with the lightest meson, the pion, in which case our general formula is exact. 
The amplitude for the whole process Map is related to the pion-a subprocess through: 

iMap =< V I -ig-r.NNl'o \ V > Z :TiManF^NN{K^) , (5) 

t — mi 



where p and p' are the initial and final proton momenta, q^nn is the proton-pion 
coupling constant, and Ft^nn{K'^) is the proton-pion form factor. We are here con- 
sidering the ttq case, and for tt± we simply need to multiply Map by an isospin factor 
of \/2. Then after squaring this result we get: 

I ^Trpp I = Q-kMnQ I F-kNN I • (6) 

Hence the pion distribution inside a proton is then: 



K 



U^.U^,K-)=x I ^JJ^^^^.\F^NN{K')\'dK^ (7) 



3:2 m2 



where a.„j\fN = q'^nnI^'^- ^^ small x, the pion distribution is then proportional to x. 
That is, it will have a limited contribution for the events of our interest. The same is 
true for the transition p —>■ ttA, and also for and other scalar or pseudoscalar meson. 

2.2 Vector meson distributions: 

Here we will study the vector meson distributions, and show that these are the most 
important contributions at the small x region of our interest. In fact, most of the 
effect comes from the omega meson (proton in the final state), and the rho (nucleon 
and delta isobars in final state). 

We will need to generalize the equivalent photon approximation to the meson 
case. In this paper we just quote the results, both for spin-1 and spin-2 particles. 
Details will be presented elsewhere 0. 

(a) Omega Emission: The omega meson has purely vector interaction, so the 
matrix element for proton to proton-omega emission is: 

iM^pp =< p' I -ig^NN lii\p> F^NN{K'^)e^- (8) 



After squaring, we get: 



where the last factor comes from the completeness relation of the omega meson's 
polarization vector, after making the equivalent meson approximation. The light-like 
vector k is given by A; = (P, 0,0,— P), where P is just a frame choice parameter, 
which will not appear in the final answer. Taking the small x and small K^ limit, we 
get: 






Introducing this result in our expression for the meson distribution we finally get: 
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U^.A-.K^) = '^ \^^M^^,dK\ (10) 



(b) Rho emission: The rho meson has a vector interaction term and a "helicity- 
fiipping" term. The matrix element for proton to proton-rho is: 

tM^ = {-tg.NN < P'b> > -^ < p'WnP > qu}F,MN{K^). (11) 

For the proton to neutron-rho case an isospin factor of \/2 should be inserted in the 
above expression. After squaring we get: 






2 _ 'i-fv yiv -r Hip) 2 J pNN-^^ n , j-i f t^2\\2 



\M,,,\' = bz^(W + -'j^)\FpNN{KT, (12) 



and therefore for the rho distribution the following result is obtained: 



(V) K^ J.1 



'2 



fp^pop - / 2^r^2 I ^\FpNN{K ) \dK 

7TX J miyK'^ + mi) 

Q P^ P' 

Q p p \ pi 

where ctpTvjv = qInn/^'^ ^"^^ '^pNN = I'pNn/^'^- Numerically it turns out that the 
helicity-fiipping contribution is much smaller than the vector part contribution. 

The proton can also emit a p~ and a A++. The p -^ p^A^ transition probability 
is related to the previous one by an isospin factor of 2/3, and the p -^ p+A° by a 
factor of 1/3. There are no c or cj meson coupling to the p ^ A transition, because 
the spin of the proton is 1/2, of the A is 3/2, while the a and uj have spin zero. 
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The NAp emission matrix element can be written as 

■ fpNA ' 



rrir 



iq^ < A" I 7^7^^ I p > e^ - zg^ < A'^ I 7^7^ \ p > e^ tpNA, (14) 



where tp^A = l^ l^ 3 is the isospin factor mentioned above. Using the completeness 
relation for the isobar particle [R[: 



J2u,,{X)u^{X) = -if' + mA) 

A "^ 



m 



2 , , , If^P'u - luP'f, 



ruA 



(15) 



and the completeness relation for the p in the equivalent meson approximation, we 
get an expression for the emission matrix element squared whose value in the small 
X limit is: 



MppA |' = 






K\K' + m')[K'' + (3mi + M')K']\FpNA\ ■ (16) 



So finally the distribution function becomes: 



K 



Jp^pAyX, K 



tpNAC^pNA f K\K^ + 3mi + Ml 



^ p 



-2m2 ,r^2 



{K"^ + ml) 



FpNA{K')\'dK' (17) 



where apNA = Ipna/^'^ ^^^ the isospin factor t^^p = 1)2/3,1/3 for p ^ p A++, 
p —>■ p°A+, p -^ p^A°, respectively. 

(c) K* Emission: Here the final state can be either a E or a A. Only the latter 
gives a contribution comparable with those of the uj and p cases. In the same way as 
before we get: 



f: 



p^K'A 



^K*NA 



K^ 



K' 



2\\2jt^2 



vrx J m\,{K'^ + m]^,] 



\Fj,,MKTdK' 



+ 



(t) K 



K' 



2\\2iT^2 



„(^) 



where ay^^ = g^NA/^^^ and a^^A = f^NA/^^^ 



4:7ix J m'j^,M^{K^ + mj^,] 



(T) 



\FK*NA{K')\'dK 



2.3 Spin 2 meson distributions 

These contributions are progressively smaller due to the higher mass of the spin 2 
particles. Thus we will only consider the lightest spin 2 particle, the /2(1270), and 
we will leave its coupling constant as an adjustable parameter. 
The f2NN matrix element can be written as: 

tMf.MN = '-^^ < p' I V^n P > 0^.F/,iViv(i^'), (19) 

where the symmetric tensor 0^j, is the spin-2 particle field. The vertex V^'^ is0: 

yt.v ^ li^gMg^ + (gt^Y + S^'Y), (20) 

with s^ = p^ + p'^. 

After squaring, averaging over initial spins and summing over delta spins, we get: 

where the h.^^p„ factor comes from the equivalent particle approximation, and is given 
by: 

_ 1 1 K^ 

J^liv.pa — AdfJ-pQvcF + gfiaQvp) + „ (OrP^)"^ 






(2xP2)4 \^2 m2 ' mj ' "''^"''^"'^"^- 



Thus, in the small K^ and small x limit we get: 

'^^=^^' - ^;^^^2 ~ ^ + ;;5\4M2 + ir2)l^/.A^Mi^ )l ' (22) 
and therefore we obtain for the /2 distribution function: 



(23) 



2.4 Parton Distributions in Mesons 

From the above calculations we have obtained various distribution functions fm/pi^i Q^] 
in the equivalent meson approximation. In order to compare our results in Eq. (1) 
with the experimental data from HERA[0], we need the parton distribution functions 
F^(/3, Q"^) for each of the mesons we have considered (vr^, 7r+, p^, p°, p+, K*^, /2). 

The parton content of mesons is presently poorly known. Thus we will assume 
here that the valence, sea and strange quark distributions of all the mesons mimic 
the parton distributions of pions. For the pion's distribution functions, we use the 
GRV parameterization [0. Then each of the FJ"(x, Q^) can be written in terms of 
three independent functions: (1) the valence quark distribution: xv'^{x,Q^), (2) the 
light sea-quark distribution: xq'^{x,Q^), and (3) the strange sea-quark distribution: 
xs''{x,Q^). That is: 

F^{x, Q") = C'^xy'^ix, Q") + C^xq{x, Q^) + C^x-s{x, Q^). (24) 

The coefficients for the charged mesons m = (7r+, p~ , p^) are (C]^, C™, C^) = (|, ^' f): 
and for the neutral mesons m = (7r°,p°, /2) are {CY',C^,C^) = (^, ^' §)• For the 
special case of the strange meson /C*+, for simplicity we will assume a parton distri- 
bution F^{x, Q"^) similar to that of the charged mesons. 

3 Results 

There are three parameters in our model, namely, the momentum squared cut-off Q^ 
(which should be around Ag^^), the tensor- meson coupling constant aj^^N = ^4^^ , 
and the form factor momentum cut-off scale Kj^nn- We have chosen the values 
i^2 = QMS[GeV% ff.NN = 4, and Aj.nn = l.29[GeV]. 

In Fig. 2 we compare the results of our model with the experimental data from 
HERAg]. We have plotted the results for F^{x, Q\ (3) for Q^ = 12, 28, 75[GeV^] and 
/3 = 0.1, 0.4, 0.9. Form the plots we can see that the x-dependence of -F^(x, Q^, [3) can 
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Figure 2: Comparison between the predictions of the peripheral meson model and 
experimental data[0. 



be explained from the x-dependence of the equivalent meson content of the proton. 

4 Conclusions 

We have presented here an equivalent-meson model in order to explain the observed 
large rapidity gap events at HERA. We have seen that the peripheral meson content 
of the proton can explain the x-dependence of the measured F^{x,Q^, P) structure 
functions. Thus, from standard nuclear physics knowledge and the parton distribu- 
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tion of mesons, the essential features of this class of events can be explained. Our 
model implies the existence of interesting final states in the forward baryon, includ- 
ing helicity- flipping proton and isosp in-changing states (e.g. A"*""*") in the forward 
direction, which should be interesting to observe in large rapidity gap events. 
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